function y=fun7_1(x1,x2,x3,x4); 
y= -1.6692+3.3383*x1+3.3373*x2+3.3400*x3+3.3387*x4;